! ***************************************
FINISH
/CLEAR
/FILN, E61
/TITLE, MECHANICAL STRESS ANALYSIS OF SKIRT SUPPORTING ZONE OF HYDROGENATION REACTOR
/UNITS,SI
! ********* 趨*********
PI=3.1415926                            ! Բ
TONGT_DI=2813                           ! Ͳھ
TONGT_T=87                              ! Ͳں
TONGQ_H=615                             ! Ͳ嵽ȹ֧дĸ߶
FENGT_T=52                              ! ͷں
FENGT_DI=2833                           ! ͷھ 
XYHR1=20                                ! ڲԲ뾶
XYHR2=20                                ! Բ뾶
H1=1500                                 ! ȹ߶
H2=1000                                 ! Ͳ峤
XQUNZ_H=70                              ! HζͼȹӲֱ߸߶
QUNZ_T=22                               ! ȹں                                 
FENGT_DO=FENGT_DI+2*FENGT_T             ! ͷ⾶
PRESI=8.83                              ! ѹ
PRESA=PRESI*(TONGT_DI**2)/((TONGT_DI+2*TONGT_T)**2-TONGT_DI**2)  
	                                ! ѹͲ˲ƽ
MASS=270000                             ! 豸
PRESG=4*MASS*9.8/(PI*((TONGT_DI+2*TONGT_T)**2-TONGT_DI**2)) 
                                        ! 豸Ͳ˲Ӧ 
EXX1=2E5                                ! ϵĵģ
NUXY1=0.3                               ! ϵĲɱ
!****************ǰ***************************
/PREP7
ET,1,PLANE82                            ! 嵥Ԫ
MP,EX,1, EXX1                           ! ϵĵģ
MP,NUXY,1, NUXY1                        ! ϵĲɱ
KEYOPT,1,3,1                            ! 趨Գѡ
!****************ģ***************************
K,1,0,0 	                        ! ĵ
K,2,TONGT_DI/2+TONGT_T,0                ! Hζͼ϶
K,3,TONGT_DI/2,0                        ! Hζͼ϶ڲ
K,5,TONGT_DI/2+TONGT_T,-TONGQ_H         ! HζͼȹӲ
K,6,TONGT_DI/2+TONGT_T-QUNZ_T,-TONGQ_H	! HζͼȹӲڲ
K,7,KX(6)-XYHR1,-(TONGQ_H-XQUNZ_H)	! HζͼȹӲԲ
L,2,3                                   ! 
L,2,5
L,5,6

CIRCLE,7,XYHR1                          ! ȹӲԲԲ
LDELE,5,7,1,1                           ! ɾԲ

L,6,4                                   ! 
K,7,0,KY(8)                             ! ¶ؼ7ΪԲǸߵĵ
L,7,8                                   ! HζͼԲߵĸ
CIRCLE,1,FENGT_DO/2                     ! ڶӦԲ
LDELE,7,9,1,1                           ! ɾԲ

LSBL,10,6                               ! ߶
LDELE,7,,,1                             ! ɾ߶

L,8,10                                  ! ٴԲߵĸ
K,30,KX(10)-XYHR2,KY(8)-XYHR2           ! 帨ؼ
K,40,KX(8),KY(30)                       ! 帨ؼ
 	
L,40,30                                 ! ɸȷڲԲ
CIRCLE,1,FENGT_DO/2+XYHR2               ! ԲȷڲԲ


LSBL,7,12                               ! 󽻵ȷڲԲ
LDELE,13,14                             ! ɾ߶
KDELE,30                                ! ɾĹؼ
LDELE,9,11,1,1                          ! ɾ߶

CIRCLE,14,XYHR2                         ! ڲԲԲ
LDELE,7,,,1                             ! ɾĹؼ߶
LDELE,10,11,,1                          ! ɾĹؼ߶

L,1,14                                  ! ɸ
L,9,14                                  ! ɸ
K,5000,KX(14),KY(14)                    ! ٶڲԲ
LSBL,9,7                                ! ߶
LDELE,12,,,1                            ! ɾ߶

L,7,14                                  ! ɸ
LSBL,6,10                               ! ߶
LDELE,12,,,1                            ! ɾ߶
LSBL,8,7                                ! ߶
LDELE,6,,,1                             ! ɾ߶

K,1,0,0                                 ! ¶
CIRCLE,1,FENGT_DI/2                     ! ڱڵԲ
LDELE,6,8,1,1                           ! ɾ߶

K,40,KX(3),KY(7)                        ! 帨ؼ
L,3,40                                  ! ɸ
L,12,14                                 ! ɸ
LSBL,12,6                               ! ߶
LDELE,8,,,1                             ! ɾ߶

LDELE,3                                 ! ɾ߶
K,100,KX(5),KY(5)-H1                    ! ȹ߶εĹؼ
K,200,KX(6),KY(6)-H1                    ! ȹ߶εĹؼ
L,100,200                               ! ȹ߶
K,300,KX(3),KY(3)+H2                    ! Ͳ߶εĹؼ
K,400,KX(2),KY(2)+H2                    ! Ͳ߶εĹؼ
LDELE,1                                 ! ɾ߶
L,300,400                               ! Ͳ߶ 
L,300,11                                ! Ͳڱ߶
L,400,2                                 ! Ͳ߶
L,5,100                                 ! ȹ߶
L,6,200                                 ! ȹڱ߶
K,800,KX(2),KY(11)                      ! Ӧ·ؼ
K,900,KX(5),KY(4)                       ! Ӧ·ؼ
K,6000,KX(300),KY(300)-400              ! Ӧ·ؼ
K,7000,KX(400),KY(6000)                 ! Ӧ·ؼ
L,1,5000                                ! ɸ
LSBL,13,15                              ! ߶
LDELE,2                                 ! ɾ߶   
A,12,14,10,11,300,400,2,800,900,5,100,200,6,4,8,9,7     ! 

!***************************************
AESIZE,1,10                             ! 嵥Ԫߴ
MSHKEY,0                                ! ʷ
AMESH,1                                 ! ʷ
LREFINE,4,,,1                            ! Բ
LREFINE,11,,,1                           ! ڲԲ
FINI 
!*******************************************
/SOLU
SFL,6,PRES,PRESI                        ! ڱʩѹ
SFL,16,PRES,PRESI                       ! ڱʩѹ
SFL,17,PRES,PRESI                       ! ڱʩѹ
SFL,1,PRES,PRESG-PRESA                  ! Ͳ˲ʩƽغ
DL,3,1,UY,0                             ! ȹ׶߶ʩλԼ
DL,7,1,UX,0                             ! ͷԳʩXλԼ
SOLVE                                   !          
FINISH



!*******************************************
/POST1
PLNSOL,S,INT,0,1                        ! Ӧǿͼ
!****************Ӧ·1--1****************
NODE11=NODE(KX(800),KY(800),KZ(800))    ! ڽڵ
NODE12=NODE(KX(11),KY(11),KZ(11))
PATH,SET1,2                             ! ·Ƽ·       
PPATH,1,NODE11                          ! ݽڵ㶨·
PPATH,2,NODE12                          ! ݽڵ㶨·
PRSECT,-1                               ! г·Ի
!****************Ӧ·2--2****************
NODE21=NODE(KX(4),KY(4),KZ(4))
NODE22=NODE(KX(900),KY(900),KZ(900))
PATH,SET2,2        
PPATH,1,NODE21
PPATH,2,NODE22
PRSECT,-1
!****************Ӧ·3--3****************
NODE31=NODE(KX(10),KY(10),KZ(10))
NODE32=NODE(KX(7),KY(7),KZ(7))
PATH,SET3,2        
PPATH,1,NODE31
PPATH,2,NODE32  
PRSECT,FENGT_DI/2
FINI